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Abstract 

In this article, we propose a new approach to calculate the Nambu-Gorkov 
propagator intuitively with some linear algebra techniques in presence of the 
scalar diquark condensates. With the help of energy projective operators, 
we can obtain relatively simple expressions for the quark propagators, which 
greatly facilitate the calculations in solving the Schwinger-Dyson equation to 
obtain the gap parameters. 

PACS : 12.38.-t, 12.38.Aw, 26.60.+C 

Key Words: Diquark condensate, color superconductivity 

1 Introduction 

At sufficiently low temperature and high density, the hadrons may be crushed into 
quark matter, which may exist in the central region of compact stars, where the 
densities are above the nuclear density and the temperatures are of the order of tens 
of keV P3I21- The color superconducting quark matter may also be created in future 
low-energy heavy ion colliders. The essence of color superconductivity is the pairing 
between the quarks when there exists an attractive interaction at the Fermi surface, 
just like the BCS mechanism in QED. Perturbatively, the quark-quark interaction 
is dominated by the one-gluon exchange, which is simply the one-photon exchange 
in QED multiplying the group factor T^jT^\ the group factor can be decomposed 
into an anti-triplet 3 with attraction and a sextet 6 with repulsion, 

T^jT^l = - ^^^^ {SuSjk - SijSki) + ^ {SiiSjk + Sij5ki). 

The quarks, unlike the electrons, have color and flavor as well as spin degrees of 
freedom, many different patterns of pairing are possible, which lead to the copious 
phase structures in color superconductivity theory, such as the color-flavor locking 
(CFL) state |3j, gapless color- flavor locking (gCFL) state 2SC state, gapless 2SC 
(gS2C) state 2SCus state and Larkin-Ovchinnikov-Fudde-Ferrell (LOFF) state 
|6j, etc. Intense theoretical works have been done to investigate the different re- 
gions of the QCD phase diagram in past years The scaling formula relating the 

"'^Corresponding author; E-mail, wangzgyiti@yahoo.coni.cn. 



transition temperature and the energy gap to the chemical potential and the run- 
ning coupling constant have been derived systematically at sufficiently high density, 
which may serve a rigorous proof of the existence of the color superconducting phase. 
At the region of moderate baryon density or below, it is too difficult to warrant a 
first principle investigation with the present analytical and numerical techniques. 
Furthermore, the charge neutrality condition as well as the (3 equilibrium have to be 
taken into account to form bulk matter inside the compact stars, which induce mis- 
matches between the Fermi surfaces of the pairing quarks and puts severe constraint 
in determining the true ground states. We have to resort to the effective actions, 
for example, the QCD-inspired Nambu-Jona-Lasinio (NJL) like models jH], or solv- 
ing the Schwinger-Dyson equation (SDE) for the gaps consistently In solving 
the SDE, the manipulations in color-flavor-spin spaces is tedious and cumbersome. 
In this article, we propose a new approach to calculate the quark propagator in 
presence of the scalar diquark condensates. 

The article is arranged as follows: we introduce the Nambu-Jona-Lasinio model 
for illustrating the formation of the scalar diquark condensate in section II; in section 
III, we propose a new approach to calculate the Nambu-Gorkov propagator explicitly 
in presence of the scalar diquark condensates, which is of special importance in 
solving the Schwinger-Dyson equation; section IV is reserved for conclusion. 

2 Nambu-Jona-Lasinio model for diquark conden- 
sate 

At low temperature and high ( or moderate) baryon density, the quarks can pair 
with each other at their Fermi surfaces with sufficiently strong attractions and form 
Cooper pairs. The colored Cooper pairs ( represented by the field A") have inter- 
actions with the quarks, its vacuum expectation values ( the diquark condensates ) 
induce the energy gaps for the quarks in the superconducting phases. The formation 
of the diquark condensates can not be derived from the first principle of QCD, we 
resort to the QCD-inspired NJL model or extended-NJL model to illustrate this sub- 
ject 0. The NJL model or extended-NJL model have the same symmetries as QCD 
and can describe the spontaneous breakdown of chiral symmetry in the vacuum and 
its restoration at high temperature and high density. 

For simplicity, we can consider the two-flavor Lagrangian density with SU{2)l x 
SU{2)ji chiral symmetry, the extension to the three-flavor NJL model with SU (3) l x 
SU{3) f( chiral symmetry is straightforward and easy, 

C = q{t^^d^ - mo + 7V)g + Gsliqqf + {l^l^rgf] + GD[{tfee'j,q){zqee\,q'')], (1) 

here q'" = Cif", (f = q'^C with C = ^7^7°. The rriQ is the current quark mass and 
the ^ is the chemical potential. The quark field q = q^^i with a = 1,2,3 are color 
index and i = 1,2 are flavor index, r = (r^,T^,r^) are Pauli matrixes, {eY^ = e'^^ 
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and (e**)"^ = e"^'' are totally antisymmetric tensors in the flavor and color spaces 
respectively. 

We can introduce the auxiliary fields and perform the standard bosonization to 
obtain the following linearized model in the mean-field approximation, 

L = q{iYd^ -rno + ^^n)q - q{a + i^^TTi)q 

-l^*\ifee'^,q) - ^^^^iiqee'^.q"^) + • • • , (2) 

with the bosonic fields 

A*" ~ i(f ee^^^q , A*'' ~ iqee'^^^q^ , cr qq , tt ~ iq^y^rq . (3) 

The conditions {a) ^ and indicate that the chiral symmetry and the 

color symmetry are spontaneously broken respectively. Here only the red and green 
quarks participate in the condensate, while the blue quarks (with color index 6) do 
not. We can introduce the constituent quark mass 

m = mo + ((T). (4) 

The spontaneously chiral symmetry breaking induced by the formation of chiral 
condensate in the vacuum and the dynamical chiral symmetry breaking induced 
by the strong color attraction at low energy region can all result in constituent 
masses for the quarks, for the u and d quarks, m„^(i ~ 330MeV, for the s quark, 
ms ~ 500MeV. The constituent quark mass m changes according to the variations 
of the baryon density, in this article, we take the mass matrix for the m, d and s 
quarks to be 

m \ 

m , (5) 

m + j 

here the m denotes the degenerate mass and the the mass breaking term bm can 
be taken into account by modifying the corresponding chemical potential for the s 
quark. 

3 Nambu-Gorkov propagator with diquark con- 
densates 

We introduce the Nambu-Gorkov formulation for the basis q^^ qgu, qrd, Qgd, 

T T / C C C C \ 

^ — [Qruj Qguj Qrd, Qgd'i Qrm Qgm Ird^ Qgd) t 

* = {Qru,qgu,qrd,qgd;Q^u^qgu,Q?d^Qgd) , (6) 
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and write down the inverse Nambu-Gorkov propagator in momentum space, 



G" 
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[Go]-' 



(7) 
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(8) 



(9) 

(10) 



At moderate baryon density, the quark masses have contribution from the chiral 
condensate and can not be neglected. We have to use the following energy projective 
operators for massive free particles ||9j, which are used extensively in solving the 
Bethe-Salpeter equation (TU], 



m] 



(12) 
(13) 



with the energy E = + m? and the basic properties. 



A±(p)A±(p-) 
A±(p)A^(p-) 
A+(p) + A_(p) 
7oA±(p)7o 



a±(p1 , 

, 

1 , 

a±(p1 . 



With those projective operators, we re-write the matrix elements for the normal 
quark propagator in the following form which can greatly facilitate the calculations 
as the terms concerning the operators A+ and A_ decouple from each other when 
inverting the matrix. 



±1-1 



Gq 



7o(po - 
7oAh 



ET)A+ + 7o(po + ^^)A- 
7oA- 



+ 



(14) 
(15) 
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here = E ± /j,. 

In the following, we invert the matrix in Eq.(7) with some linear algebra tech- 
niques to derive the Nambu-Gorkov propagator G intuitively. For the CFL and 

gCFL states, the is a 72 x 72 matrix in the color-fiavor-spin spaces, it is very 
difficult to obtain the Nambu-Gorkov propagator G with simple expressions, the 
present approach can greatly facilitate the calculations. 
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(16) 



1 - 
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Gn A+ 



[GoV 



G+A-Go-A+ 



GoA+G+A- 



1 - 



G+A- 



GoA^ 







[G+]-'-A-Go-A+ 

[Go]"'-A+G+A- 

G+ -G+A-G- \ _ f G+ E- 
G- j ~ I S+ G- 



-Gq A+G+ 



(17) 



(18) 
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with 

(7+ = ^ , (19) 

G- = ^ . (20) 

In the 2SC and g2SC phases with the basis "if in the color-flavor spaces (see Eq.(6)), 
the matrix A"*" and A~ (see Eqs.(8-10)) are skew diagonal, the following matrix is 
strictly diagonal, 



G+A-G'oA+ 

GqA+G+A- 



(21) 



the diagonal and shew-diagonal matrix elements decouple from each other in Eq.(16). 
2 We can read out the matrix elements Gfu , G^^, G% ,G%, , Sj^,.^ 

, 'S^u,rd from Eqs. (17-20) straightforwardly with the help of the energy projective 
operators. For example, 

G: 



'ru 



[Ga;:-A(-^75) [Gj]^,A(-^75) 

1 



7o(Po - i??JA+ + 7o(Po + EtJ^_ + A75(-^ + 



7oA+ + 7 ;;tt7 — r^ToA- 



{po + E^J (po - E%) - A2 ' ^ (po - EfJ {po + EJ,) - A 



= GUru^gd) 



gd 

--ru,gd - ~ 777—^ T^X^ 777 + 7 / 7:^7 777^- 



iPo + EtJ {po - E%) - A2 ^ (po - Ef^) {po + EJ,) - A^ 
"± ( 

' 'ru,gd ' 'ru,gd\ 



- (ru^gd). (22) 



We can write down other matrix elements analogously with simple substitution. 

In the CFL and gCFL phases, the scalar diquark condensates can be parameter- 
ized as 

{q?C-f,q^) ~ A,e"^%,, + A,E^^^e,,, + A,e^P%,, , (23) 



^The results in Eqs. (18-20) certainly can be derived by the standard linear algebra manipula- 
tions, in this article, we propose to separate the diagonal and shew-diagonal matrix elements in 
Eq.(16) intuitively, the matrix in the denominator in Eq.(17) is block diagonal, for the 2SC and 
g2SC states, the blocks G',]"A~G|^A+ and GqA^GqA~ are strictly diagonal, the matrix can be 
inverted easily. The general strategy be, choose the suitable basis for the matrixes A+ and A~ to 
diagonalize or skew-diagonalize them in the color-flavor spaces, the resulting blocks Gq Gq A'^ 
and GqA^GqA~ are strictly diagonal or mostly diagonal with some non-diagonal blocks. The 
matrix G~^ can also be diagonalized by analytical plus numerical calculations, which is performed 
in Ref.|ll|. Furthermore, there are also other works using the energy projective operators to invert 
the propagator G~^, for example, Ref.|12|. 
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the condensates are Lorentz scalars, antisymmetric in Dirac indices, antisymmetric 
in color (with the strongest attraction between quarks) and antisymmetric in flavor. 
The gap parameters describe d — s, s — u and u — d Cooper pairs, respectively. We 
can introduce the Nambu-Gorkov formulation in the three-flavor case, 

,c „c „c ^ ^c-^ ^ 

bu' Qgs-i Qbd) ■ 

(24) 

The inverse Nambu-Gorkov propagator is a 18 x 18 matrix in the color-flavor spaces, 

\ 



^ — i^Qru-i Qgd-i Qbs-i Qrdi Qgu, Qrs, Qbu, Qgs, Qbd, Qr 
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(25) 



(26) 



(27) 



v 

where A is 3 x 3 block matrix, B, C and D are 2x2 block matrixes. The matrix 



" Go A+G+A- 



G+A- 



(28) 



is block diagonal, the matrix GjA~Go A"*" and Gq A+GjA" are diagonal for the six 

basis {rd, gu,rs,bu, gs,bd) and have non-diagonal blocks concerning the matrix A 
for the three basis {ru,gd, bs). The diagonal matrix elements of the inverse matrix 
of 



■G0-A+ 







[Go] 



A 



-G+A- 



(29) 
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are the inverse of the corresponding matrix elements. We can read out the matrix 
elements G^^ ^G^u^Gf^^G^^^Gf^^ G^^ directly from the Eqs. (17-20) with the help of 
the energy projective operators. For example, 

1 



G. 



K];;-A3(-i75)[G^],„A3(-z75) 



r7oA+ + 7 ^^v^^^^;T^-r^7oA- 



{Po + E%) (po - E^^) - Ai + (po - E^,) (po + E^:) - A 



Gtu = G^^{gu^rd), 



^± ^ -^^375 ~ -iA375 ~ 

(Po + E^,) (po - ) - Ai + ^ (po - E^,) (po + - Ai - ' 

^%,rd = S^^ ff„(r(i ^ £?li) . (30) 
We can write down the corresponding ones for the Gf^, G^, Gf^, G^^, (,^, 

^%M^ -bd,9s analogously. 

For the non-diagonal blocks concerning the A in Eqs. (28-29), in general one 
may expect a linear transformation for the basis {ru, gd, bs) can result in diagonal 
matrixes for both the block A and the corresponding [Gq"] however, it is not 
the case. The color superconducting phases may exist in the core of compact stars, 
where the bulk matter should satisfy the charge neutrality condition and the /3- 
equilibrium. The gravitation force is much weaker than the electromagnetic and the 
color forces, any electric charges or color charges can forbid the formation of bulk 
matter. For a neutral two-flavor color superconductor, the quark chemical potentials 
in color and flavor spaces can be expressed in terms of baryon chemical potential 
electrical chemical potential fig, and color chemical potentials /^s and 

2 

while for the three-flavor quark system, the chemical potentials are taken as 

2 8 5m^ 

From above equations for the chemical potentials, we can see that they are far from 
being equal. The charge neutrality condition puts strong constraints on both the 
two-flavor and three-flavor systems and results in the g2SC and gCFL phases. In the 
gCFL phase, the up block A for the matrix A^ in Eq.(25) and the corresponding one 
for [Gq ] ^ in Eq.(27) can not be diagonalizcd with the same basis. Fortunately, we 
can solve the SDE for the gaps Ai, A2, A3 with the six basis {rd, gu, rs, bu, gs, bd). 

The explicit and simple expressions for the quark propagators are of great im- 
portance in solving the gap equations and determining the gap parameters A, Ai, 
A2, A3. 

G-\p) - G^\p) = G{k) rf D^/^{k - p) , (33) 



/^m,/3j = fJ'SapSij — fieQijSa/S + l^zT^p^ij H -=ll^T^i^5ij ^(^i3^j3<^a/3 ■ (32) 
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where 











(34) 



The gluon propagator has complex energy and chemical potentials depen- 

dence, include the effects of Landau damping and Debye screening induced by the 
medium Furthermore, the contributions come from the scalar diquark conden- 
sates should also be taken into account in the color superconducting phases. The 
gluon self-energy has been investigated in the g2SC and gCFL phases. The results 
show that the gauge bosons connected with the broken generators have imaginary 
Meissner screening masses ^2], the chromomagnetic instability should be taken into 
account in a consistent analysis. The explicit expressions for the quark propaga- 
tors can greatly facilitate the calculations, especially in performing the Matsubara 
frequency summation. At low baryon density, the coupling constant can not be 
factorized out, g'^D^^ should be approximated by some phenomenologically poten- 
tials just as the conventional SDE at zero baryon density ^M,- 

In the CFL state with the same chemical potential (or in the limit /i oo) and 
degenerate mass m {6m ^ 0), we can perform a linear transformation to diagonalize 
the matrix A and diag{GQ^^, Gq^^, Gq^^) simultaneously. 
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With the new basis (xi, X2, x^), the matrix A takes the following form 
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are strictly diagonal for the nine basis 
(xi, X2, xs, rd, gu, rs, bu, gs, bd), we can read out the G^ for the basis {xi, X2, x^) from 
Eqs.( 17-20) directly. 
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The block matrix for the G and S with the basis (r-u, gd, hs) can be written as 



G± I O/^i /^i /^i \ 

G^-Gj G^ + 2Gt G^-Gt , (39) 
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(40) 



From above matrixes, we can see that for the normal propagators, due the presence 
of the scalar diquark condensates, all the matrix elements (qaiQ/sj) are non-zero, while 
only the diagonal elements in the color-flavor spaces are non-zero in the ordinary 
vacuum. 



4 Conclusion 

In this article, we propose a new approach to obtain the Nambu-Gorkov propagator 
intuitively with some linear algebra techniques in presence of the scalar diquark con- 
densates. With the help of energy projective operators, we obtain relatively simple 
expressions for the quark propagators, which can greatly facilitate the calculations 
in solving the SDE to obtain the gap parameters. 
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